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Abstract. We propose simple model of colloidal sessile drop desiccation. The model describes 
correctly both evolution of the phase boundary between sol and gel inside such a drop and the final 
shape of the dried film (deposit). The model is based on mass conservation and natural assumption 
that deposit (gel phase) prevents flows and evaporation. 
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1. Motivation 

During last decade, the drying of particle-laden sessile drops attracts much attention from the sci- 
entific community because of numerous applications. First of all one should refer to the important 
surface-patterning technique known as evaporation-induced self-assembly (EISA) or evaporation 
driven self-assembly (EDSA), ink-jet printing, coating technologies, spotting technologies for bio- 
assays, etc. 

As a rule, the particle-laden sessile drops desiccate with pinned contact line. The outward flow 
occurs inside the drop and carries the particles to the drop edge. The ring-like deposit remains on 
the substrate ^M- 

There may be two very different situations when a particle-laden sessile drop desiccates. First, 
the particles inside such a drop may interact with each other only pure mechanically (impacts). 
In this case, the deposit forms a porous medium. This medium prevents neither bulk flow inside 
it nor evaporation from its free surface. Second, the particles may form chemical bonds. In the 
last case, a phase transition may occur {e.g., sol to gel). The new phase blocks any flows, the 
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solvent evaporation from such a phase is negligible. The proposed models mainly deal with the 
first situation E SI 13 CI S HX51 H31 H21 UHl ESI EH 133 - Only several models treat the deposit as 
impenetrable for flows and inhibiting evaporation lQ3lfl4ll28ll2"9ll . Nevertheless, the simulation of 
desiccated colloidal drops with phase transition is extremely important for high-throughput drag 
screening ES, biostabilization OH, and medical tests flMEDEl. The models lfl3l fl4ll28l l29ll 
utilize the sets of rather complicated PDE's. We propose alternative very simple but efficient 
model. 

There are three different modes of colloidal sessile drop desiccation IfToll : 

1 . t G t D , where t G is the gelation time and t D is the desiccation time. The gelled phase occurs 
near the drop edge and moves inward while the central area of the drop remains liquid. 

2. t G pa t D . The gelled skin covers the free drop surfaces. This thin shell cannot prevent 
evaporation of the solvent. The buckling instability occurs. 

3. t G <C t D . The phase transition from sol to gel in the whole bulk of the drop is almost 
instantaneous. The gelled drop loses solvent via evaporation very slow. 

The mode of our interest is the first one. In this particular case, the desiccation process can be 
divided into several stages. 

1. Initial single-phase liquid stage (Figure [DO- The whole drop is a sol. The outward flow 
carries suspended particles to the drop edge until volume fraction of the suspended particles, 
$, reaches critical value, $ g . Note that particle-enriched region is extremely narrow whereas 
the particle volume fraction in central area of the drop is almost constant along its radius E51 

2. Intermediate two-phase stage (Figure EEb,c). Gelled ring near the drop edge is formed and 
growths to the drop center. The volume fraction of the colloidal particles is constant both 
inside the gelled foot, $ g , and sol, except rather narrow area near the phase front Il28ll2~9ll . 

3. Final single-phase solid stage (Figure [TJl). The gelled deposit loses the rest of moisture 
and cracks. The real fluids of interests (e.g., biological fluids) may content both suspended 
particles and dissolved substances. In this case, the dendritic crystals may occur in the central 
area of a sample [fTll24|. 

In the present article, we focus our efforts on the second stage only. The construction of this 
paper is as follows. In Section [2]1 we describe our model. In Section [3] the obtained results are 
discussed. Finally, we summarize our results and conclude this paper. 
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2.1. Main assumptions 

We consider a colloidal (nano- or microparticle-laden) sessile drop deposited on horizontal solid 
impermeable substrate. The drop is axisymmetric. In our consideration, we suppose 

1. The liquid is incompressible. 

2. The contact line is pinned during the desiccation. 

3. Free surface of the liquid phase is a spherical cap. This assumption is valid for the sufficient 
small drops when Bond number is much smaller than unity and rather slow flows when 
capillary number is small. 

4. Evaporation from the gelled phase is negligible and doesn't effect on hydrodynamics inside 
liquid phase. This assumption is supported by the experimental data |fT2~ll2~3l . 

5. Any hydrodynamical flows inside gelled phase are absent. This statement can be concluded 
from the simple observation that gelled ring cracks (see, e.g., JT), i.e. there is not any flow to 
prevent moisture loss even by very slow evaporation. 

6. The drop apex height as well the drop mass decrease almost linearly up to full gelation (T). 

7. The phase boundary is vertical. Volume fraction of the colloidal particles is in the 
liquid (sol) phase and constant $ g in the solid (gel) phase. This assumption is supported by 
our simulations [|28l . 

8. Sedimentation is negligible. 

2.2. Governing equations 

In the cylindrical coordinates, free surface of the liquid phase, z(r, t), is described by the equation 



where hf(t) is the height of the phase boundary between sol and gel, rj(t) is the radial coordinate 
of the phase front, H(r, t) is the profile of the free surface above liquid (sol) phase and H (t) = 
H(0,t) (see Figure Sh). 

The conservation of solute yields 



2. Model description 
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All quantities, $(£), rj(t), hf(t), are. the functions of time, t. We omit arguments here and below 
to simplify notation. 

The conservation of solvent yields 



dt r^ /+ 6 + 2 r f 



) +r f ^(2h f (* a -*) + H (l-$))+ 



°(i/ 2 + rJ) + (l-$)r^ = -- y J( r| t)Wl+^-J rdr, (2.2) 



dt 



where p is the density of solution. 

Sum of (12.11) and (12.21) gives the equation for the solution 





Note that right hand side of (12.21) . (12.31) is proportional to drop mass variation in time 



dm f , . / / <9/A 2 , 

- = -2„j J(r, t )Jl + (-\r d r. 



We denote 



? J J M)y X +(§) rdr = ^). 



In spite of the fact that during evaporation, drop mass decrease looks almost linearly IfTl |3~1 l23l 
l32l . deviation of the evaporation rate from constant is fairly visible [23]. We suppose that evapo- 
ration rate is proportional to the area of free liquid surface and decreases when volume fraction of 
colloidal particles increases 

A(t) = J rj(t) (l - ^) • (2-4) 

To our best knowledge, the experiments regarding density of vapor flux above two-phase desiccated 
sessile drops are not performed yet, so reasonableness of our choice (12.41) can be examined only 
by simulations. 

Mass conservation for the solution can be rewritten as 

TT drt ldH 0/rr r, 9 . 9 dhf . ,„ „. 

r > H °nr + 2nr (H ° +r < )+r 'nr = - A (2 - 5) 



We suppose 

d_ 

dt 



(h, + H ) = V (l - ^ , (2.6) 
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where constant V is the rate of apex height decrease for a drop of pure solvent. Our assumption is 
supported by experiments (see, e.g., lfTll23lO. 

The set of equations (12.11 ), (12.51 ), (12.61 ) is underdetermined. Additional assumptions are re- 
quired. 

From physical point of view, contact angle between sol, gel, and air, 6, have to be constant for 
arbitrary position of the phase boundary, if the volume fraction of colloidal particles inside liquid 
phase is fixed. Obviously, for the high concentrated colloids ($ ~ $ 9 ), 6 have to tend to zero. 
Contrary, for the dilute colloids (<3> <C O — > O , where O is the contact angle between gel, 
air, and pure solvent. We assume 

9(f) = Go (l - ^) • (2.7) 

In fact, this speculation is valid for the static equilibrium of the phases. If the phase boundary 
moves, not only surface tension, but viscous forces should be involved into consideration. Un- 
fortunately, to our best knowledge, there is not any published measurements of O , thus, it is an 
adjustable parameter in our model. 

At the point A with the coordinates (rf,hf), angle 7 between tangent line to the liquid phase 
and horizon, angle f3 between tangent line to the gelled phase and horizon (see Figure [2]) can be 
found from the relation 



2 




substrate 



Figure 2: 1 and 2 are the tangent lines in point B to the free surface of gelled phase and liquid (sol) 
phase, respectively, a is the angle between these lines, f3 and 7 are the angles between tangent 
lines and horizontal direction. 

It is clear, that 

. . dhf 

tan/3(t) = -— - 1 (2.9) 
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and 6 = 7T — a = n — (ix — 7 + /?) = 7 — /3 (see Figured]). Hence, from (12.81 ) and (12.91 ) yields 

dhf ( nt \ 2r f H (t) \ dr f ( ^. . 2r f H (t) \ „ n 

-d y + Un6it) r}-H°\o,t) ) = -i ( tane(i) - ) • (2M) 

where 0(£) is defined by Equation l2.7l 

Set of equations (12.11) . (12.51) . (12.61) . (12.71) . (12.101) is complete and allows calculate the shape of 
two-phase colloidal drop and volume fraction of colloidal particles inside sol at any time. 

It is convenient to measure all lengths in units of initial drop radius (radius of its constant base). 
In this case < r < 1. 



2.3. Parameters of the model 

It is rather difficult to extract necessary information from the published experimental data. 

To estimate parameters of the model, we use, particularly, the data about desiccation of drops 
of a phosphate buffered saline solution (pH = 7.4, ionic strength / = 0.2 M) of bovine serum 
albumin (BSA) fl}. For the concentration of BSA 40 gl _1 , initial drop mass is about 15 mg, initial 
height is about 1.5 mm, and drop radius is about 2.5 mm. Thus, initial ratio of height to radius is 
about 0.6. Central zone of the deposit is flat and its height is 36 ± 2 /mi. Thus, ratio of final deposit 
height to the initial drop height is about 0.024. For the concentration of BSA 60 gl~\ central part 
of the deposit is about 60 /im, while the rim is about 80 /im. 

For the drop with radius 1.75 mm and concentration of albumen 100 gl _1 , the heights of the 
dip and rim are about 15 and 35 fim, respectively 0. 

For the drops of native blood plasma and radius 3-3.5 mm, the heights of the dip and rim are 
about 10 and 20-40 yum, respectively |fTT|. 

In our computation, initial volume fraction of colloidal particles <3> varies from 0.04 to 0.4, 
and ratio of initial drop height, H(0, 0) to its radius varies from 0.06 to 0.6. 



3. Results 

Figure [3] demonstrates our simulations for H(0,0)/R = 0.06, $(0) = 0.1. Evaporation rate is 
Jo = 2.2 ■ 10 -4 , velocity of apex decreasing is Vb = 4.22 • 10~ 4 , and adjustable parameter is 
®° = ifo • Different profiles correspond (from up to down) to the initial time (t = 0), t = to/4, 
t = to/%t = 3£d/4, and t = to- Final profile of the gelled film is close to the experimentally 
observed 0][EI ED • 

Dynamics of volume fraction of the colloidal particles looks like reported by Sobac and Brutin 
in 11231 . The experiment |[23l shows rather sharp decreasing of the drying rate, \dm/dt\, when the 
system changes from two-phase (sol and gel) to one-phase (gel). Unfortunately, our computations 
of drying rate don't demonstrate the same excellent agreement with the experiment, i.e. drying 
rate decreases monotonically. We suppose the results of simulations may be improved by means 
of more accurate definition of (I2.4I ). 
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Figure |4] shows evolution of the drop during desiccation if initial volume fraction of the sus- 
pended colloidal particles is 0.04 and the drop is rather thick (H (0, 0)/R = 0.6). These parameters 
correspond to the experiments by Annarelli et al. (T). Adjustable parameter 6 is n/8. Instants 
are the same as in Figure |3] 




r 



Figure 4: Dynamics of shape. Thick drop. 



4. Conclusion 

Proposed semiempirical model give simple explanation of the colloidal sessile drop desiccation. 
The results obtained for the parameters corresponding blood serum are consistent with experi- 
ments. The model allows to describe correctly the final shape of the dried film (deposit). 
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